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d Abstract 



Plasmodium of Physarum polycephalum is a single cell visible by unaided eye. On a 
non-nutrient substrate the plasmodium propagates as a traveling localization, as a 
compact wave-fragment of protoplasm. The plasmodium-localization travels in its 
originally predetermined direction for a substantial period of time even when no 



p , gradient of chemo-attractants is present. We utilize this property of Physarum lo- 



calizations to design a two-input two-output Boolean logic gates {x, y) — )■ {xy, x + y) 
and (x, y) — )■ {x, xy). We verify the designs in laboratory experiments and computer 
simulations. We cascade the logical gates into one-bit half-adder and simulate its 
functionality. 

Keywords: Physarum polycephalum, logical gate, adder, unconventional computer, 
chemical computer, biological computer 
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'^ 1 Introduction 



/\ Plasmodium of Physarum polycephalum is a single cell with many diploid nu- 

5^ clei. The cell is visible by naked eye and can grow up to meters when properly 

cared for. The plasmodium feeds on microscopic food particles, including mi- 
crobial life forms. The plasmodium placed in an environment with distributed 
nutrients develops a network of protoplasmic tubes spanning the nutrients' 
sources. 

In its foraging behavior the plasmodium approximates shortest path [13], com- 
putes planar proximity graphs [5] and plane tessellations [16], exhibits prim- 
itive memory [15], realizes basic logical computing [17], and controls robot 
navigation [18]. The plasmodium can be considered as a general-purpose com- 
puter because the plasmodium simulates Kolmogorov-Uspenskii machine — 
the storage modification machine operating on a labeled set of graph nodes [4] . 
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In 2004 Tsuda, Aono and Gunji [17] demonstrated in laboratory experiments 
realisation of Boolean logic negation and conjuction by plasmodium of Physarum 
polycephalum. In 2004 Adamatzky and De Lacy Costello esablished in numer- 
ical simulation [3] and chemical laboratory experiments [9] that by colliding 
localized excitations, or wave- fragments, in excitable chemical medium one can 
implement functionally complete set of logical gates. We merge approaches [17] 
and [3,9] in present paper. We adapt concepts of collision-based computing [2] 
to realms of Physarum behaviour, and develop experimental prototypes of 
two-input two-output Boolean logical gates. 

The paper is structured as follows. Methods of cultivating and experimenting 
with Plasmodium of Physarum polycephalum are described in Sect. 2. In Sect. 3 
we provide experimental evidence of 'ballistic' behavior of traveling Plasmod- 
ium localizations. Experimental Physarum gates are discussed in Sect. 4. In 
Sect. 5 experimental results are supported by numerical simulation of propa- 
gating localizations. The gates are cascaded in one-bit half-adder in Sect. 6. 
Importance of non-nutrient substrate for gate implementation is highlighted 
in Sect. 7. 



2 Materials and Methods 



Plasmodium of Physarum polycephalum is cultivated in large plastic boxes, 
on a wet paper towels and fed with oat flakes. Experiments are conducted in 
round Petri dishes (9 cm in diameter) and rectangular Petri dishes (12 cm x 
12 cm). Channels and junctions physically representing logical gates are cut of 
a non-nutrient 2% agar plates (Select agar, Sigma Aldrich). The dishes are kept 
in room temperature (c. 25") in darkness. Images of plasmodium propagating 
in Petri dishes are taken by Epson Perfection 4490 scanner, resolution 600. 
Colors are enhanced by increasing saturation and contrast. 

We use two-variable Oregonator model to numerically simulate propagation 
of Plasmodium localizations. Our choice and details of the model are outlined 
below. 

Localized excitations in sub-excitable Belousov-Zhabotinsky (BZ) medium be- 
have similarly to pseudopodia of P. polycephalum [6,7]. Sources of nutrients are 
chemo-attractants for plasmodium, gradients of shade are 'photo-attractants' 
for excitation waves in BZ medium. In [7] we shown how to navigate trav- 
eling localizations and growing parts of plasmodium by spatial configuration 
of attractants. We adopt the analogy developed in [7] and simulate propa- 
gating Plasmodium using two-variable Oregonator equation [10] adapted to a 
light-sensitive BZ reaction with applied illumination [8]: 



ot e u + q 

dv 

In framework of BZ reaction the variables u and v represent local concen- 
trations of activator, or excitatory component, and inhibitor, or refractory 
component. With regards to plasmodium of P. polycephalum activator, u, is 
analogous to concentration, or 'thickness', of the Plasmodium's cytoplasm at 
the propagating pseudopodium. Inhibitor, v, combines several factors, when 
Plasmodium is concerned. These factors include rate of nutrients consumption, 
byproducts of biochemical chains ignited by signals on photo- and chemo- 
receptors, and concentration of metabolites released by the plasmodium into 
its substrate. 

Parameter e sets up a ratio of time scale of variables u and f , g is a scaling 
parameter depending on rates of activation/propagation and inhibition, / is 
a stoichiometric coefficient. Constant is a rate of inhibitor production. In 
light-sensitive BZ represents rate of inhibitor production proportional to 
intensity of illumination. In terms of plasmodium represents rate of inhibitor 
proportional to concentration of nutrients, metabolites, illumination, chemical 
repellents. See detailed comparison of BZ and Physarum in [6,7]. 

To integrate the system we use Euler method with five-node Laplacian oper- 
ator, time step At = 0.001 and grid point spacing Ax = 0.25 (equivalent to 
0.6 mm of physical space), e = 0.0243, / = 1.4, q = 0.002. The parameter 
characterizes excitability of the simulated medium: the medium is excitable 
and exhibits 'classical' target waves when = 0.05 and the medium is sub- 
excitable with propagating localizations, or wave-fragments, when = 0.0766. 



3 Ballistic of Physarum localizations 



Proposition 1 Given cross-junction of agar channels and plasmodium inoc- 
ulated in one of the channels, the plasmodium propagates straight through the 
junction. 

We experimentally found that that plasmodium propagates under its own mo- 
mentum when no gradients of repellents or attractants are applied. In example 
shown in Fig. la the plasmodium is inoculated in northmost part of the north- 
south channel. The plasmodium has the only option — to propagate south — 
because there is no gel substrate further north. Thus an internal 'momentum' 
is formed. No food sources are applied on the substrate to attract the plas- 
modium. The Plasmodium propagates straightforwardly by itself. It does not 
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Fig. 1. Experimental examples of plasmodia moving under their own momenta: 
(a) no sources of chemo-attractants, (be) oat flakes are placed in south and east (b), 
and south and west (c) ends of channels. 

branch at the junction and moves till it reaches south end of the north-south 
channel (Fig. la). 



In 21 out of 28 trials the Plasmodium exhibits clear 'ballistic' behavior and 
propagates straight through the junction as under its own momentum. In 7 
out of 28 trials the plasmodium turns into other channels or branches into 
several channels at once. Adding sources of attractants — assuming they are 
balanced over the channels (i.e. if there is an attractant in east-west channel 
there should be one in north-south channel) slightly speeds up propagation of 
Plasmodium and does not change overall statistics of the plasmodium propa- 




q p 

Fig. 2. Geometrical structure of Physarum gates Pi (a) and P2 (b): x and y are 
inputs, p and q are outputs. 

gation. Examples are shown in Fig. Ib-d. Thus in experiments described in 
Sect. 4 we always used oat flakes to stimulate the Plasmodium growth. 



4 Physarum gates 



Geometrical structure of gates Pi and P2 is shown in Fig. 2. We experimented 
with various shapes of agar and found that the most suitable templates are 
those shown in Fig. 2. Input variables are x and y and outputs are p and q. 
Presence of a plasmodium in a given channel indicates Truth and absence 
— False. 

Each gate implements a transformation from {x^y) — )■ {p,q)- Experimental 
examples of the transformations are shown in Fig. 3. Plasmodium inoculated in 
input y propagates along the channel yq and appears in the output q (Fig. 3a). 
Plasmodium inoculated in input x propagates till junction of x and y, 'collides' 
to the impassable edge of channel yq and appears in output q (Fig. 3b). 

When Plasmodia are inoculated in both inputs x and y they appear in both 
outputs p and q (Fig. 3c-e). In some cases plasmodia originated in different 
inputs avoid each other (Fig. 3cd) and thus head towards different outputs. 
In other cases the plasmodia merge in a single plasmodium but nevertheless 
branch towards different outputs (Fig. 3e). There are no strict rules on re- 
pelling and merging and often initial repelling between two plasmodia can be 
followed by merging (Fig. 3fg). 



Typically for biological substrates, plasmodium of Physarum polycephalum 




Fig. 3. Experimental examples of transformation (x, y) — >• (p, q) implemented by 
Physarum gate Pi. (a) (0,1) ^ (0,1), (b) (1,0) ^ (0,1), (cHe) (1,1) ^ (1,1). 
(f)-(g) snapshots of Pi gate taken in 12 h interval. 

is sensitive to experimental conditions. The Plasmodium sometimes deviates 
from general scenario when implementing transformation (x,y) — )■ {p,q)- Re- 
sults of experiments with Physarum gate Pi are shown in Fig. 4a. 

Plasmodia do not cancel each other at once. Therefore if at least one of the 
inputs is '1' we expect to see '1' at one of the outputs. Input scenario (1, 1) 
is straightforward: in six out of seven experiments plasmodia appear on both 
outputs. Thus we obtain transformation (1, 1) — )■ (1, 1) (Fig. 4a). 

Plasmodium inoculated in input y (while input x is empty) will appear in 
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Fig. 4. Experimental data on transformations (x, y) — )■ {p, q) implemented by 
Physarum gates (a) Pi and (b) P2- Values x = 1 and y = 1 are represented by 
Plasmodia inoculated in inputs x and y, respectively. Values p = 1 and q = 1 are 
represented by plasmodia reaching outputs p and q, respectively. Frequency of each 
particular scenario (a, b) — )■ (c, d) is presented by fraction: denominator is a total 
number of experiments for {x, y) = {a, b) and numerator is a number of experiments 
completed with output tuple {x,y) = {c,d). 

output q in 17 out of 22 experiments. Thus the transformation (0, 1) — )■ (0, 1) 
is realized by plasmodium in over 70%. 

Input combination (1,0) gives us less stable results: in nine our of 13 ex- 
periments (69%) the Plasmodium reaches output q. The plasmodium enters 
output p in three of 13 experiments, and the plasmodium branches inside both 
outputs in one experiment. Nevertheless the transformation (1, 0) — > (0, 1) is 
realized by plasmodium in well over half of experiments (Fig. 4a). 

Finding 1 Plasmodium of Physarum polycephalum implements two-input 
two-output Boolean gate {x, y) — )■ {xy^ x + y) with reliability exceeding 69%. 

Experimental snapshots of plasmodium propagating in the gate P2 are shown 
in Fig. 5. Taken input x is empty, plasmodium placed in input y usually (see 
statistics in Fig. 4) propagates directly towards output q (Fig. Sab). Plasmod- 
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Fig. 5. Experimental examples of transformation (x, y) — )• (p, q) implemented by 
Physarum gate P2: (a) (0,1) — )• (1,0), plasmodium is inoculated in input y, no 
oat flakes present; (b) (0,1) — )• (1,0), oat flakes are placed in both outputs; 
(c) (1,0) — )• (0,1), no oat flakes present; (d)-(e) two snapshots of transformation 
(1,1) —7- (0,1), taken with 11 h interval, oat flakes are placed in both outputs; 
(f) (1,1) — )• (0,1), oat flakes are placed in both outputs; (g)-(h) transformations 
(1, 1) — )• (0, 1) are less pronounced than in previous examples, however we see that 
output q is more extensively occupied by plasmodium than output p. 



ium inoculated in input x (when input y is empty) travels directly towards 
output p (Fig. 5c). Thus transformations (0,1) — )■ (1,0) and (1,0) — )■ (0,1) 
are implemented. 





(a) (b) (c) 

Fig. 6. Time lapsed images of localized excitations, wave- fragments, traveling in 
channels of gate Pi filled with excitable medium. Dynamics of the excitable medium 
gate Pi is shown during implementation of transformations (a) (0,1) — )■ (0,1), 
(b) (1, 0) -^ (0, 1), (c) (1, 1) -^ (1, 1). The transformation (1, 1) -^ (1, 1) is simulated 
for initial excitations in channels x and y positioned at equal distance from their 
meeting point. 

The gate's structure is asymmetric, a:-channel is shorter than y-channeL There- 
fore the Plasmodium placed in input x usually passes the junction by the time 
Plasmodium originated in input y arrives at the junction (Fig. 5d-f). The 
y-plasmodium merges with x-plasmodium and they both propagate towards 
output q. Extension of gel substrate after output q does usually facilitate im- 
plementation of the transformation (1, 1) — )■ (0, 1) (Fig. 5d-h). 

Frequencies of various input-output transformations occurred in experiments 
are shown in Fig. 4b. Plasmodium inoculated in input y will reach only output 
p in 21 out of 29 experiments. Transformation (1, 0) — ?■ (0, 1) takes place in 
16 out of 27 experiments. Transformation (1,1) — )■ (0,1) occurs in 13 out of 
21 experiments. 

Finding 2 Plasmodium of Physarum polycephalum im,plem,ents two-input 
two-output gate {x, y) — ?■ (x, xy) with reliability exceeding 59%. 



5 Simulation of Physarum gates 



Experimental Findings 1 and 2 are confirmed in numerical simulation using 
two- variable Oregonator model of a sub-excitable medium (see Sect. 2). To 
represent value '1' in input channel x or channel y we generate an excitation 
near entrance of the channel. Two wave-fragments are formed. One travels 
outside the gate, another travels towards outputs. 



Wave-fragments in sub-excitable media are notably unstable, they do keep 
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Fig. 7. Time lapsed images of localized excitations, wave- fragments, traveling in 
channels of gate P2 filled with excitable medium. Dynamics of the excitable medium 
gate P2 is shown during implementation of transformations (a) (0,1) — )■ (1,0), 
(b) (1,0)^(0,1), (c) (1,1)^(0,1). 

their shape only for a short period of time. Then the fragments either shrink 
and annihilate or expand unlimitedly. During simulation of gates Pi and P2 
we manually adjusted parameter e (see Sect. 2) to keep wave-fragments from 
collapsing and expanding. Figures 11 and 12 illustrate dynamics of e during 
simulations of gates Pi (Fig. 6) and P2 (Fig. 7). 

Let us look at the time lapsed images of wave-fragments propagating in gate 
Pi (Fig. 6) and gate P2 (Fig. 7). Scenarios where only one input is excited 
are straightforward. When y = 1 in gate Pi is initialized the wave-fragment 
propagates southward along the channel yq (Fig. 6a). Wave- fragment initiated 
in input x {x = 1) oi gate Pi propagates along input channel x and collides 
to the boundary of channel yq (Fig. 6b). The wave-fragment recovers after 
collision to the boundary, travels along the channel yq and appears in the 
output q (Fig. 6b). Wave-fragments behave similarly in situations of input 
tuples (1, 0) and (0, 1) in gate P2. They propagate straight along their original 
input channel and reach outputs opposite to their entry points (Fig. 7ab). 

In input scenarios (0, 1) and (1,0) size of propagating wave-fragment was not 
enough for the fragment to branch into output channel p of gate Pi (Fig. 6ab). 
When two wave-fragments are initiated, x = 1 and y = 1, they collide at the 
junction of input channels x and y. The wave-fragments merge into a single 
larger- wave fragment. This new fragment propagates towards output q and 
also expands into output channel p (Fig. 6c). 

In gate P2, inputs x = 1 and y = 1, wave-fragment originated in x- input 
arrives at the junction between channel xp and yq a bit earlier then wave- 
fragment originated in y-input. Therefore y-wave collides into refractory tail 
of a;- wave. The y-wave annihilates (Fig. 7c). Such development is phenomeno- 
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x^ 



xy + xy 

Fig. 8. Scheme of one-bit half-adder made of gates Pi and P2- Inputs are indicated 
by arrows. Outputs xy + xy and xy are sum and carry values computed by the 
adder. Outputs and x + y are byproducts. 

logically identical to Physarum gate P2 (Fig. 5de) with the only difference that 
Plasmodium 'wave' originated in input y is not annihilated but merges with 
Plasmodium originated in input x. 



6 Simulated one-bit half-adder 



One-bit half-adder is a logical circuit which takes two inputs x and y and 
produces two outputs: sum xy + xy and carry xy. To construct a one-bit 
half-adder with Physarum gates we need two copies of gate Pi (Fig. 2a) and 
two copies of gate P2 (Fig. 2b). Cascading the gates into the adder is shown 
in Fig. 8. Signals x and y are inputted in P2 gates. Outputs of P2 gates are 
connected to inputs of Pi gates. 

We did not manage to realize one-bit half-adder in experiments with living 
Plasmodium because the plasmodium behaved differently in the assembly of 
the gates than in isolated gates. Simulation of the adder using Oregonator 
model was successful (Fig. 9). To simulate inputs x = and y = 1 we initiate 
wave-fragments at the beginning of channels, marked y and arrow in Fig. 8. 
The wave-fragments propagate along their channels. The waves do not branch 
at the junctions with other channels because we keep the wave-fragments 
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Fig. 9. Time lapsed images of localized excitations, wave- fragments, traveling in 
channels of one-bit half-adder filled with excitable medium. Dynamics of excitations 
is shown for input values (a) x = and y = 1, (b) x = 1 and y = 0, (c) x = 1 and 
y = l. 
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localized by varying parameter e (Fig. 9a). 

For input values x = 1 and y = wave-fragments are originated in sites 
marked x and arrow in Fig. 8. The wave-fragment started in left x-input 
channel propagates towards x + y-output of the adder. The wave-fragment 
originated in right x- input channel travels towards xy + xy (Fig. 9b). 

When both inputs are activated, x = 1 and y = 1, wave- fragment originated in 
left y-input channel is blocked by refractory tail of wave-fragment originated 
in left x-input channels. The wave-fragment traveling in right x-input channel 
is blocked by tail of wave-fragment traveling in right y-input channel. The 
wave-fragments representing x = 1 and y = 1 enter top-right gate Pi and 
emerge at its outputs xy and x + y (Fig. 9b). Thus functionality of the circuit 
Fig. 9 is demonstrated. 



7 Discussion 



We established experimentally and in numerical simulations that plasmodium 
of Physarum polycephalum realizes basic logical operations on a geometrically- 
constrained non- nutrient substrate. We designed two types of Boolean logic 
gates, both gates have two inputs and two outputs. The gates implement 
transformations {x,y) — )■ {xy,x + y) and {x,y) — )■ {x,xy). We shown how the 
Physarum gates can be assembled into a one-bit half-adder. Functionality of 
the adder is illustrated using two-variable Oregonator model. 

Our designs are based on the interactions between traveling localizations: 
Plasmodium localizations propagating on a non-nutrient substrate and wave- 
fragments propagating in a sub-excitable medium. Similarities between the 
Plasmodium localizations and wave-fragments are discussed in details [6]. We 
stress that things go absolutely differently on a nutrient-rich substrate (corn 
meal agar) and fully excitable chemical medium. Plasmodium inoculated in 
any point of the nutrient-agar gel gate propagates in all channels (Fig. lOa-c). 
An excitation wave initiated at any point of excitable medium gate spreads 
all over the gate (Fig. lOd). We can conclude therefore that it is impossible to 
implement logical functions with plasmodium of Physarum polycephalum on 
a nutrient-rich substrate. 

Reliability of experimental Physarum gate is quite low: 69% for gate Pi and 
59% for gate P2. This is because behavior of plasmodium is determined by 
too many environmental factors — thickness of substrate, humidity, diffusion 
of chemo-attractants in the substrate and in the surrounding air volume, and 
physiological state of plasmodium during each particular experiment. Increas- 
ing reliability of Physarum gates might be a scope of further studies. 
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(a) 



(b) 




Fig. 10. Snapshots of plasmodium propagating in gate Pi on 2% corn meal agar 
(a~c) and numerical simulation (d): (a) inputs x = and y == 1, see the gate's 
scheme in Fig. 2, (b) inputs x = 1 and y = I, (c) inputs x = 1 and y = 0, (d) time 
lapsed images of excitation wave propagating in gate Pi for inputs x = 1 and y = 0, 
simulation of case (c). 
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Appendix A 



Figures 11 and 12 sliow liow parameter e (see description of tlie model in 
Sect. 2) is changed during simulation of gates Pi and P2 in a sub-excitable 
medium. We also show dynamics of the medium's activity. At each step of 
simulation we calculate activity a* = |L|~^ ■ J^x^l'^x ^^ ^ sum. of values u for 
each node of the grid L normalized by total number of nodes |L|. 
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(c) 

Fig. 11. Dynamics of parameter e, dotted line, and activity a, solid line, during 
operation cycle of gate Pi for input tuples (a) (0,1) (see Fig. 6a), (b) (1,0) (see 
Fig. 6b), (c) (1,1) (see Fig. 6c). 
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(c) 



Fig. 12. Dynamics of parameter e, dotted line, and activity a, solid line, during 
operation cycle of gate P2 for input tuples (a) (0,1) (see Fig. 7a), (b) (1,0) (see 
Fig. 7b), (c) (1,1) (see Fig. 7c). 



